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Introduction 

Characterizing entanglement is one of the challenging issues that has been fostered by 
the development of quantum computation in the past few years. Beyond the obvious 
interest for the foundations of quantum mechanics, the motivations increased as the 
role played by entanglement in the power of quantum algorithms was made clearer. 

As random states are entangled with high probability, they play an important 
role in the field of quantum communication, and appear in many algorithms, such 
as quantum data hiding protocols [Ij or superdense coding [2]. Various quantum 
algorithms were proposed to generate random states: algorithms based on the 
entangling power of quantum maps have been proposed in [3]; chaotic maps were 
considered for instance in [H [5] , pseudo-integrable maps in [6] . It was also proposed 
to generate random states by construction of pseudo-random operators [7j, or by 
sequences of two-qubit gates [8]. The efficiency of all these algorithms is measured 
by their ability to reproduce the entangling properties of random states. Thus, it is 
crucial to have some measure of entanglement. 

The problem of quantifying entanglement is a difficult issue, and a number of 
entanglement measures have been proposed (see e.g. the review [9]). However for 
pure quantum states, bipartite entanglement, that is the entanglement of a subset 
of the qubits with the complementary subset, is essentially measured by the von 
Neumann entropy of the reduced density matrix [10] . More convenient to use are 
the linear entropy or the purity, which are linearized versions of the von Neumann 
entropy. Purity or linear entropy were used to estimate the entangling properties of 
chaotic quantum maps (for instance in the Baker's map [5j or in kicked tops ^11)). 
or entanglement growth under random unitary evolution [12]. It has been used as a 
reference to compare the accuracy of various entanglement measures |13j , or to measure 
dynamical generation of entanglement in coupled bipartite systems |14l . Recently, it 
was shown that purity for a pure quantum state could be expressed as a function of 
the inverse participation ratio, thus enabling to connect entanglement (as measured 
by the purity) to localization properties of quantum states |151 116j . 

In most of these works, the entanglement properties of the systems that are 
studied are compared to those of random pure states. This is usually done by resorting 
to numerical computations, based e.g. on Hurwitz parametrization 17j of random 
states. Of course, a comparison directly based not on numerics but on analytical 
formulae describing the entanglement properties of random pure states would be most 
desirable. In this paper we derive such formulae for the purity and the linear entropy. 

Random pure states can be realized as column vectors of random unitary matrices 
drawn from an ensemble with unitarily invariant Haar measure |18j . Equivalently they 
can be realized as vectors with coefficients given by independent random complex 
Gaussian variables, rescaled to have a norm equal to 1 (see e.g. [19]). Various aspects 
of the entanglement properties of random pure states have been studied in previous 
works: the distribution of G-concurrcnce has been calculated in [20]; the average von 
Neumann entropy had been obtained in 21j; the moments for the purity distribution 
in random states were calculated analytically in j22j , as well as approximate moments 
of Meyer- Wallach entanglement (a multipartite entanglement measure based on the 
purity). The aim of this paper is to provide analytic expressions for the probability 
density distribution of the purity (or the linear entropy) for bipartite random pure 
states. Explicit analytical formulae are derived for the smallest bipartitions of the 
Hilbert space into a p-dimensional and a g-dimensional spaces: p = q ^ 3 'm Section[3l 
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p = 3 and any q in Section 31 p = g = 4 in Section [51 A different method allowing to 
obtain formulae for p = A and g > 4 is explained in Section O The method is general, 
and similar calculations would yield expressions for higher values of p and q, although 
it is not sure whether these formulae would take any nice and compact form. 



1. Probability density distribution of the purity. 

Let us consider a state belonging to the Hilbert space ® C for some integers 
p, q with p < q. Suppose 5* admits the following Schmidt decomposition |23| 

p 

i=l 

where \ai) and \bi), I < i < p, are respectively orthonormal bases for and C. The 
purity R of the state can be written in terms of Schmidt coefficients as 

= (2) 

1=1 

The linear entropy is expressed in terms of the purity by the simple relation 

^LW = ^(l-i?(*)). (3) 

p-1 

For random pure states obtained as column vector of random matrices distributed 
according to the unitarily invariant Haar measure (CUE matrices), the Schmidt 
coefficients are characterized by the following joint distribution [19]: 

Pix^,...,xp)=A n n 4"''Mi-E^H (4) 

l<i<j<p 1<2<P \ «=1 / 

for Xi G [0, 1]; (5 is the Dirac delta function, and A is the normalisation factor 

Uo<,<p-iii-j-mp-jy' ^ ' 

The purity distribution function is then given by 

/ JL. \ ( JL. \ 

, (6) 



P{R)^A jdP^Vi^f Jl ^-^'(^ [ l-^x, j (5 [i?-^ 

•^0 l<i<P V i=l / V 4=1 



where we have introduced the Vandermonde determinant V^(x) — Y\i<ji^i ~ ^j) 
with X = {xi,...,Xp). The distribution P{Sl) of the linear entropy can be 
straightforwardly deduced from P{R) using Eq. ([3]): 

PiS,) = ^-lp (l- P-ls,) . (7) 



P 

As the purity vanishes outside the interval [1/p, 1], the distribution of linear entropy 
is supported by [0, 1]. 

In the next sections we evaluate ^ in the cases p = 2, 3, 4. We first note that 
the expression ([6]) can be simplified to 



P{R)=Ap\ f\pxVix) Yl x'^+''-^S ll 5 Ir~J2 

•'^ l<i<p \ 1=1 / V i=l 

where d = q— p, using a transformation detailed in [22] . 
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2. Distribution of the purity for a 2 x gr bipartite random state. 

In the case p = 2,q > 2, the analytic expression for the probabihty distribution P{R) 
can easily be obtained analytically directly by integration of (l6|). It reads 

(2g-l)! 



^^^^ 2«-i('Z-l)!(9-2)! 
for 1/2 < i? < 1, and otherwise. 



(1 - i?)9-V2i?- 1 



(9) 



3. Distribution of the purity for a 3 x 3 bipartite random state. 

In the case p = 3 and g = 3, Eq. ([6]) reads 



PiR) 



4 JO 



d^:>cV{x.)x2xlS 1 




(10) 



The domain of integration is the intersection of the cube [0, 1]'^, the plane X1+X2+X3 = 
1 and the sphere of radius VR centered on the origin. This intersection is a circle if 
i < i? < i (see Fig. [T]left), or sections of a circle if ^ < R < 1. We first make the 



change of variables x 
/ - 



V 



2 

ye 

V6 



ox, where O is the orthogonal matrix 




1 

V2 



t , 



(11) 



In the new coordinates (Xi,X2,X3), the plane has equation X3 = l/-\/3 and the 
sphere X\ + X\ + X| — R. We first integrate over X3 and make a change of variables 
to polar coordinates {Xi — pcos9,X2 — psin9). Let us set = if | < R < \, 
(f) — arccos(l/\/6i? — 2) if 5 < i? < 1, and r = \/R — 1/3. The domain of integration 
in the plane — 1/ ^/i is represented in Fig. [T] right for two different values of R. 
After integrating over p we obtain 





(12) 



f ) 1 sin(3^). 



The domain of integration for 9 is I?^ 



u 



[0, 27r] \ Cfjj, where 
27r , 27r" 
3 3 



u 



47r 



47r 

T 



(13) 



The integrand in Eq. (|T2|) can be expanded as a sum of cos k9 and sin A:6'. Terms of 
the form sin k9 yield zero upon integration; terms of the form cos k9 yield 



d9 cos k9 



d9 [1 + con - 



2kTT 



Akn 



cos k9. 



(14) 



which is zero unless 3|A;. Nonzero contributions therefore necessarily come from the 
cos [9 + cos^ (0 + ^) sin(3^) term in (fT2)l . Keeping only terms of the form cos 3kt 
in the expansion of this term we get 

PiR) = 70V3r^ (27r - dcj) + sin(60)) . (15) 
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Replacing cfi and R by their value we finally obtain 



P{R) = 70V3 



2n [R- 



— arccos 



(16) 



-(i?- |)V6i?-3, 5<i?<l. 



It is interesting to check whether this formula allows to recover the moments derived 
in [22]. The value of (i?") is given (see [22) by 

pl{pq — 1)! nl 



E 



(pQ + 2n — 1)\ ^ nilno^- ■ 



i=l ■ l<i<i<p 

(correcting by a factor p\ the Eq. (10) of where this term had been erroneously 
forgotten). The first moments of P{R) as calculated from Eq. ITBl) yield values that 
are precisely equal to those obtained from Eq. ([T7)) . 





Figure 1. Domain of integration for p = 3. Representation in for ^ < < | 
(left), and in the plane xi + X2 + X3 = 1 (right, thick lines). 



4. Distribution of the purity for a 3 x g bipartite random state, 9 > 3. 



In the case of a random vector belonging to a Hilbert space of dimensions p x q 
with p = 3 and 9 > 3, the term V{xi,X2,X3)x2x'^ in (jlOp is multiplied by a factor 
(xiX2a;3)'^, where d = q — p. The changes of variables of Section [3] yield (still setting 
r — \/ R — 1/3) an integrand that can be expanded as 



,,3 a 

V ^ fe=0 



1 

27 




COS 9 



sin 36* cos'' 36*. 




(18) 



Again in psp only terms of the form cos{3j6) contribute to the final result, and nonzero 
contributions therefore necessarily come from the 



fk {e) = cos ( + y j cos^ + Y ) sin(36l) cos'' 36* 



(19) 
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term. For a given k, the expansion of cos'^ 39 yields 



k/2-j 



cos 6j9 k even 



cos*= 361 = . Y \ (20) 

Eflo'^^' ((fc _ i)/2 _ cos(6j + 3)0 k odd, 

where Sj is the Kronecker delta. For a fixed integer t the only terms contributing 
to P(i?) in the expansion of cos (^6 + ^) cos'^ [9 + ^) sm{39) cos 3t9 sum up to 
a/3 (cos(3i -6)9-2 cos3t9 + cos(3i + 6)6*) /32. The integral of fki9) over thus 
yields 

3^3 ^ / 1 \ / k 



E ( 1 - ¥n (k/L J ^ ^^^^ (21) 



2fc+3 Zl^ I 2 V \k/2-j 

3\/3 ' f k 

2fc+3 2^ 



((fc_l)/2-j)^^+i/2(<^)' ^odd, (22) 



where we have defined, for j > 0, the function 

sin(6j-6)0 sin6j(j> sin(6j + 6)0 
^^(^)= 6,-6 6j + 6 (2^) 

(for j = and j = 1 it is understated that the limit j — > or j ^ 1 is taken). 
Summing all contributions together we finally obtain the exact expression 

(3(7-1)! M / -1 \''V5-9-R\'^"''Vd 1^^^ 



lk/2\ 



1 



^ V 2 



with again </> = for 1/3 < i? < 1/2 and arccos(l/V6it; - 2) for 1/2 < i? < 1. Here [a;J 
is the integer part of x, fc fc mod 2 and 5 is the Kronecker delta symbol. Note that the 
term Xj+ki4') can be expressed as a function of R in terms of Chebychev polynomials 
of the second kind Un, using the relation sinnfl = J7„_i(cos0) sin6'. Again, on can 
check that the expression p4)) allows to recover the moments (fT7|) . 



5. Distribution of the purity for a 4 x 4 bipartite random state. 

The treatment of the p = A case is quite similar to the previous one but calculations 
(and results) quickly get very heavy. Here we derive an explicit expression for the 
case p — q — A. The first steps of Section [3] are easily generalized (they can in fact 
be generalized in an obvious way to arbitrary p, q) . The domain of integration in ([6]) 
is the intersection of a plane and a hypersphere in R^, restricted to [0, 1]^. First we 
make a change of variables x = OX, where O is the orthogonal matrix 
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so that the plane has equation X4 = 1/2 and the sphere + X| + X| + X| — R. 
After integration over X4, we have to integrate 

f{Xi,X2,X3) = X2x\x\V {xx,X2,Xz,X/^ (26) 



over the portion of the sphere of radius r = \/ R~ 1/4 that is comprised inside 
the thetrahedron of vertices Vi = (-\/3/2, 0, 0), V2 = (1/VT2, -2/V6, 0), = 
(1/712, l/%/6,-l/V2) and V4 (l/\/T^, 1/V6, 1/V2), centered on (0,0,0). In the 
case 1/4 < R < 1/3, the sphere is entirely inside the tetrahedron: after making a 
change of variables to spherical coordinates {Xi = p cos 61, X2 = p sin 9i cos 62, X3 = 
p sin 01 sin02), the integration is trivial. Taking into account the factors 1/2 and l/2r 
coming from the integration of 6(2X4 — 1) and S{R — 1/4 — p^) the integration of (pS)) 
yields 



rfX/(X) = 47rr^-^/45045, (27) 

with X — {Xi,X2,X3). When 1/3 < i? < 1/2, the integrand can first be simplified 
by use of the symmetries of the domain of integration, which consists of the sphere 
minus the four spherical caps emerging from the tetrahedron. We note Ci , 1 < i < 4 
the cap opposite to vertex Vi. Let Ri{(p) be the rotation of axis OXi, i = 1,2, 3, and 
angle (p. The caps C2, C3, C4 can be obtained from Ci by rotations: C2 is the image 
of Ci by = i?i(— 7r/3)i?3(arccos(— 1/3)), and C3 and C4 are respectively the images 
of C2 by T = i?i(27r/3) and = i?i(47r/3). Therefore 



/ dX/(X) = f dX (/(X) + /(5X) + f{TSX) + fiT^SX)) 



(28) 



In spherical coordinates, the top cap (cap Ci) has equation p = r,0 < 9i < 
arccos , < 02 < 27r (we recall that r — \/ R— 1/4) . Performing the integrals 
over 62, and then 61, in ((28|) yields 



^X/-(x) = 9677r 571^ ^2 ^ 5057r 

804925734912 V3 9459597312 429981696^3'' 
2297r R 71" s IIti" m 19 87r ,0 ,„„, 

7-0 -\ T-*' 7^10 ^12. ^ ^13 |-2g\ 

20901888\/3 20480\/3 124416^3 20736^3 45045 
Finally, when 1/2 < i? < 1, the domain of integration can be further reduced to 
62 € [0,7r/3] by applying T and and using the symmetry 02 — > —02- The domain 
of integration is 02 G [0, 7r/3] and 0i g [(/s, 7r[ where verifies 

rcos(^ + V3/2 

= ^ . , (30) 

l\J Ar sm<y9 

and the integrand is a sum of terms of the form cos^'^' 0i sin0i cos^'" 02 and 
cos^*^"'""'^ 01 cos^'^ +^ 02. The calculation gets tedious and the steps leading to PiK) 
in this case are given in the Appendix. 

Taking into account the normalization factor Ap\ — 15!/(3!2!)^ we finally get 
P(R) = for i? ^ [1/4,1] and 



P(i?) = i^(4i?-l)-, i?e 

16 

P(i?) = -Q2(i?)-^^(4i?-l)^, i?e 
3 16 



1 1 

4' 3 
1 1 
3' 2 



(31) 
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P{R) = V6i?- 3 Qi(i?) 
4725 



— — arccos 



1 



16 



-(4i?-l)- 



arccos 



V6i? 
R 



Q2{R) 



m-i 



,Re 



1 



where we have mtroduced the following two polynomials Qi{x) — j^^^i^ ^ 
a;) (1657 + 277731a; - 2190321a;2 + 6208416a:3 - 7386066a;'* + 2913408a;5) and Q2{x) = 
^^i^(-159241 + 2178306X - 11709126a;2 + 30254796a;3 - 34540506x4 + 4864860a;^ + 

14594580a;^). Again it can be checked analytically that the moments obtained from 
the above formula for P{R) agree with those obtained from Eq. ([T7]). 



6. Probability density distribution of the purity for a 4 x q bipartite 
random pure state, q > 4. 

Similar calculations can be done in the same way for g > 4. However as the number of 
terms in the integrand gets large, it is more efficient to proceed in the following way. 
The only difference between the g = 4 and the q > 4 cases is that the integrand is 
multiplied by {xiX2X3X4)'^ , where d = q — p. It is then easy to see from Section [5] that 
the integrals appearing in the calculation for 1/4 < i? < 1/3 and 1/3 < i? < 1/2 still 
yield polynomials in r. For 1/2 < R < 1, one can show that after reducing the domain 
of integration by symmetries as in the previous section, the integrand is again a sum 
of terms of the form cos^*^ 6i sin 6*1 cos^*^ 02 and cos^'"'"'"^ di cos^*^ +* 6*2- The calculation 
therefore yields an expression similar to Eq. p5p . but with polynomials in r of higher 
degree. Therefore one can look for a P{R) (expressed in terms of r = \JR— 1/4) of 
the form Az{r) for 1/4 < i? < 1/3, AqIj) for 1/3 < i? < 1/2, and of the form of 
Eq. (j35|l for 1/2 < i? < 1, where Ai all are polynomials in r. The maximal order of 
the Ai corresponds to the order of the integrand {xiX2X^xaYx2x\x\V{xi,X2tX^,xa), 
that is (taking into account the from the Jacobian and the 1 /r of the integration of 
one of the delta functions) li + Ad. The coefficients of these polynomials are unknown 
variables. The knowledge of all moments (|17[) allows us to write down as many linear 
equations as there are unknown variables, that is 6(14+4(i). The problem now reduces 
to solving a linear system Mx — h. The vector h = {R^), {R^)i . . .} is obtained 

from (|17p . and M is a matrix of size 6(14 + Ad) whose coefficients are given by terms 
of the form 

/, r''g{r)dr, fc = 2,3,4, (32) 

where g{r) is one of the functions 1, arctany^^^, arctan \J 121^^-?, '-'^ \/4r2 — 1. 
These integrals can be calculated analytically. The solution x — M^^b of the system 
yields the coefRcients of the polynomials Ai and thus an analytic formula for P{R)- 

A similar method would yield expressions for higher Hilbert space dimensions. 
However, given what Eq. ([3T|l looks like, it is not clear whether formulae for higher 
dimensions could be cast into a tractable form. 
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Appendix 

The derivation of P{R) for p = 4, g = 4 and | < i? < 1 starts from the integrand 
obtained after having reduced the domain of integration by symmetries, and integrated 
over p and 9i S [v,7i'[. From Eq. (|30p and using the fact that 02 G [0, 7r/3[ we obtain 
an expression for sin if and cos ip which allows to express the integrand as a function 
of 02 only. Expanding all trigonometric terms and changing variables from cos 02 to 
X, we are left with a sum of terms of the form 



x'^'^VT 



(33) 



(1+8x2)13 

and 

X^a+l^l _ a;2 ^4^2(1 82,2) _ 3 
(1 + 8x2)13 ' 

where a is some integer and x has to be integrated between 1/2 and 1. The integrals 
corresponding to (j33[) give constants independent of r. The integrals corresponding 



(34) 



to (1331) can be evaluated by the change of variables t = y/ [x"^ — c)/(l — x"^) with 
c = (3 — 4r^)/(32r^). Upon integration, we obtain terms of the form 

I 8^^2 I 2 / 

Ai{r) + A2{r) arctan W _ i + ^3(?-) arctan J ^^^^ _ ^ + Ai{r)^/ Ar"^ - 1, (35) 

where Ai are polynomials in r. Replacing r by its value \ R — j and noting that 



12r2 - 1 



arctan ■ 



/ 


TT 


' 12r2 - 3 


~ 2 


/ 8r2 


TT 


4r2 - 1 


" 2' 



1 4r2 + l „ 

2 '^'^^^Q^ 12^2 _ 1 + aretan ^ = -, (36) 
we get the final result (|3T|) . 
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